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In the context of models where the dark energy density pr> is a random variable, 
anthropic selection effects may explain both the "old" cosmological constant problem 
and the "time coincidence". We argue that this type of solution to both cosmologi- 
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cal constant problems entails a number of definite predictions, which can be checked 
against upcoming observations. In particular, in models where the dark energy den- 
sity is a discrete variable, or where it is a continuous variable due to the potential 
energy of a single scalar field, the anthropic approach predicts that the dark energy 
equation of state is po = ~Pd with a very high accuracy. It is also predicted that 
the dark energy density is greater than the currently favored value Qrj ~ 0.7. An- 
l/") . other prediction, which may be testable with an improved understanding of galactic 

properties, is that the conditions for civilizations to emerge arise mostly in galaxies 
completing their formation at low redshift, z « 1. Finally, there is a prediction 
which may not be easy to test observationally: our part of the universe is going to 
^ ■ recollapse eventually. However, the simplest models predict that it will take more 

than a trillion years of accelerated expansion before this happens. 
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INTRODUCTION 



The "old" cosmological constant problem - why don't we see the large vacuum energy 
density p\ which is expected from particle physics? - and the "time coincidence" problem 
- why do we live at the epoch when the dark energy component p^ starts dominating? - 
may find a natural explanation in models where po is a random variable. The idea is to 
introduce a dynamical dark energy component X whose contribution px varies from place 
to place, due to processes which occurred in the early universe. Then 

Pd = Pa + p x 

will also vary from place to place, and the old cosmological constant problem takes a different 
form. The question is not why pa is much smaller than rj 4 , where rj is some high energy 
physics mass scale, such as the supersymmetry breaking scale i] ~ TeV, but why do we 
happen to live in a place where p^ is almost exactly cancelled by px- This line of enquiry is 
rather quantitative, since we can ask what is the probability for us to observe certain values 
of Pd ~ 10 _11 (eV A ) 4 , or what is the probability for the time coincidence. 

Explicit particle physics models for a variable px have been reviewed in [1] . Two examples 
which have been thoroughly discussed in the literature are a four-form field strength, which 
can vary through nucleation of membranes [2, 3] , and a scalar field with a very low mass [3, 4] . 
Assuming one such mechanism, and using a theory of initial conditions such as inflation, 
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one can calculate the "a priori" probability distribution V*{pD)dpD- This is defined as the 
fraction of co-moving volume which at some fiducial initial time (which we conventionally 
take to be the time of recombination) had the value of the dark energy density in the interval 
dpo- Inflation is also responsible for smoothing out the value of po over comoving distances 
much larger than the size of our presently observable universe. 

By itself, V* is not sufficient to calculate probabilities for our observations. Selection 
effects which bias the measurement of pn must be included, and the most important one 
in this case is anthropic [5-9]. 1 While \po\ may be very large in most places, there is 
nobody there to observe such extreme values. If p^ > 0, galaxy formation stops once the 
dark energy becomes dominant over the matter density. Some galaxies are seen at redshifts 
of order z ~ 5, but not much higher, indicating [9] that galaxies will not form in regions 
where p^ > (1 + zeg) 3 Po- Here, p is the matter density at the present time t , and zeg ~ 5 
is the redshift at the time t^G ~ (1 + z E G)~ 3 ^ 2 to when the earliest galaxies formed. Also, 
for a negative pp the universe recollapses on a time scale ~ [Gp^l -1 / 2 , where G is 
Newton's constant. This time should be larger than the earliest time tsi which is required 
for intelligence to develop [7, 13]. Thus, observers will only exist within a tiny "anthropic 
range" 

-(Gil;)' 1 < p D < (Gt 2 EG )-\ (1) 

It should be noted that, aside from the above minimal requirements, anthropic selection 
includes all other ways in which p E disfavours the existence of observers. For instance, in 
regions where pr> < 0, the matter density is larger than \pn\ throughout the cosmic evolution. 
If \pd\ is too large, all galaxies formed in that region will be very dense, and as a result, 
very inhospitable. This occurs also for a large p D > 0, since galaxies must form before pr> 
starts dominating. We shall come back to this issue in Sections III and V. 

The selection effect can be implemented quantitatively by assuming the mediocrity prin- 
ciple, according to which our civilization is typical in the ensemble of all civilizations in the 
universe. The probability to find ourselves in a region with given values of po is thus given 
by [14] 

dV{p D ) ccV*(p D )n C iv{pD)dpD- (2) 

Here, n civ (p£)) refers to the number of civilizations which will ever form per unit co-moving 
volume in regions where the dark energy density was equal to pr> at the time of recombina- 
tion. 2 

Needless to say, the determination of both factors in the r.h.s. of Eq. (2) leaves room for 
some uncertainties. However, we shall argue that there are reasons to be optimistic. If the 
distribution (2) is to explain both cosmological constant problems, then a number of rather 
generic predictions can be made, rendering these ideas very testable. 

In the next section, we review the calculation of the prior probability distribution V*(pd)- 
The anthropic factor n C i V (pD) is discussed in section III. In the same section, we argue that 
the anthropic approach can succeed only if the conditions for civilizations to evolve arise 
mostly in galaxies formed at low redshifts, z ~ 1. Anthropic predictions for the dark energy 



1 Anthropic selection effects associated with the possible variation of the amplitude of density fluctuations 
[10, 11] and of the baryon to photon ratio [10, 12] have also been discussed in the literature. 

2 As we shall argue, in models where both cosmological constant problems can be solved anthropically, po 
has not varied appreciably since the time of recombination, and therefore it can be treated as constant in 
time. 
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equation of state, for the energy density pp, and for the Hubble parameter h are discussed 
in sections IV, V. The prediction for the future of the universe in unveiled in section VI. 
Finally, our conclusions are briefly summarized and discussed in section VII. 

II. THE PRIOR DISTRIBUTION 

The first task in determining (2) is to estimate V*. The vacuum energy density is of 
order p\ ~ rj 4 > (TeV) 4 , and therefore pn must have a natural range of variation of order 
rj 4 or larger. Weinberg noted [9] that a function V^po) that varies smoothly on scales 
Pd ~ ff ■, should behave as a constant in the utterly narrower interval (1) - unless of course, 
the function would happen to have a zero or a pole in that interval (which would be an utter 
coincidence). This led him to conjecture that for values of po in the anthropic range the 
prior probability would be constant, 

^*(pd) ~ const. (3) 

Outside of this range the form of V* is irrelevant, because the factor n C i V vanishes. Weinberg's 
conjecture is subject to verification. As mentioned in the Introduction, V* is calculable, 
provided that the dynamics of px is known, and assuming an inflationary model which 
would determine its spatial distribution at the time of recombination. Analysis of explicit 
models shows that (3) is not automatically guaranteed [4], but it does seem to be satisfied 
in generic models. 

There are basically two reasons [1, 4] why a non-flat V* may result from the process of 
randomization of p^> which occurs during inflation (this randomization is due to quantum 
diffusion in the case where A is a scalar field, or to nucleation of membranes in the case when 
A is a four-form). The first reason is the differential expansion induced by the dark energy 
component. During inflation, the expansion rate is determined by H 2 = (8nG /3) (Vi n f + pn) ■ 
Although po is very small compared with the inflationary potential Vi n f, its effect may 
build up over time, in such a way that more thermalized volume is generated with high 
values of pr>. In this way, V^(pn) could be biased towards large values. Let us denote by 
r(A, H) the characteristic time needed for the dynamics of A to sample (at a fixed point in 
space) all values of pn within the anthropic range (Ap D ) anth . The differential expansion is 
characterized by the parameter 

q = (AH)t = (4nG/3)H-\Ap D ) anth r(X, H). (4) 

If q >> 1, then V* is exponentially steep in the range of interest. This case is ruled out 
by observations, because it predicts a very large p^, even after selection effects have been 
factored in. If q ~ 1, the distribution V* may have a moderate dependence on p D within 
the anthropic range. This dependence affects the position of the peak of the distribution for 
the observed values of po, Eq. (2), and hence it affects our predictions. While models of 
this sort are not ruled out, they require a very unnatural adjustment of parameters, since 
q is determined by a combination of rather different pieces of dynamics. Hence, we shall 
disregard this marginal possibility as non-generic. Finally, there is a wide class of models 
where q <C 1 is satisfied without any fine tuning [1, 4], and hence we shall take this to 
be the generic case. Numerical simulations confirm that in this case the bias effect due to 
differential expansion is insignificant [18]. 
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The second reason why V* may be non-flat is the following. Even if the differential 
expansion is negligible, and the prior distribution for X is flat, this does not automatically 
guarantee that the prior for p D will be flat, unless the relation between X and pr> is linear 
in the range of interest. Through this effect, it is possible to have a moderate variation of 
V*(pd) within the anthropic range. But again, this would require a contrived adjustment 
of parameters and we shall dismiss this case as non-generic (see also [19] for a discussion of 
this issue). 

As an example, let us consider the case where px = V(4>) is the potential energy density 
of a scalar field 0, 

Pd= p A + V{<f)). (5) 

The field must change very slowly on a cosmological time-scale, so that its potential energy 
behaves as an effective cosmological constant. This requires the slow-roll conditions [4] 

\V'\ < 10p D /m p , \V"\ < l0 2 PD /m 2 p (6) 

to be satisfied up to the present time (when p D ~ p , with p the present matter density). 
The constraint g< 1 on the differential expansion yields [4] 

V' 2 H 4 /GV 3 > 1. (7) 

During inflation, the scalar field is randomized by quantum fluctuations, and at recombina- 
tion it is distributed according to the "length" in field space, 

V*((f))d(f) oc d<j>. (8) 

Therefore 3 , 

V*(p D )dp D oc j^r^jj- (9) 

Thus, the flatness of the prior depends on how much V changes in the anthropic range. 
As we shall see, variations in this range may occur, but they do not bias the probability 
distribution for po in any significant way, unless we adjust some parameters specifically for 
this purpose. 



Consider a potential of the form 



v(4) = (io) 



where p 2 pA < 0, so that it is possible to have \po\ very small even if \p^\ is large. Eqs. (6) 
lead to the condition [4] 

H « 10- 120 m>Al" 1/2 - (11) 
Such a small mass parameter may seem unrealistic, but it can naturally arise, for instance, 
in a low energy effective theory with a suitable discrete symmetry [3] (for other proposals, 
see [1, 19, 20] and references therein). Note that (11) does not correspond to a fine tuning, 
but just to a strong suppression. The condition (7) translates into 

M > HlH- 2 ~ KT 169 m p > (12) 



3 Note that near the points where V'(<j>) — 0, we have po ~ A + B(j) 2 and V '((/>) <~ <j) ~ (pd — A) 1 / 2 , which 
is integrable. Hence, the zeroes of V'(<fi) are not a concern. 
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where H is the present Hubble rate, and in the last step we have used H ~ 10~ 7 m p , 
corresponding to a GUT scale of inflation. The conditions (11) and (12) leave very many 
orders of magnitude available for the parameter p, and so fine tuning is not necessary. From 
(5), 

Pz) = / t (0-0 o ) + ^(0-0o) 2 , (13) 

where 0q = —2p\/p 2 and k = p 2 0o- We are interested in the vicinity of p D = 0, where it is 
easy to show from (9) that [4] 

V*(p D )dp D oc [1 + 0{p D /p A )]dp D « dp D . (14) 

Since pr> <C pa in the anthropic range, the distribution is indeed flat to a very good accuracy. 
For contrast, we may consider the "washboard" potential 

Pd = p A + + M 4 sm((j)/ri), (15) 

where k was given above and M and rj are different mass scales. Let us assume that 4 

M 4 « Hlm.pT] ~ (rj/m p )p . (16) 

Then the field will typically be found away from the local minima, with a probability dis- 
tribution 

'P*(PD)dpD = lK+{M ^ w{<j)/r})l -V ( 17 ) 

Both k and M 4 /p should be much smaller than H^m p in order to satisfy the slow roll 
condition. In the case k ^> M 4 /i], the distribution (17) is still flat, as in (3). In the 
opposite case, where M 4 /p ^> k, the a priori distribution can have a sizable variation within 
the anthropically allowed range. If 77 <§; m p , this range is very wide in the field space, 
5(j) > Po/ K ^ m p- This means that the oscillations in V* will average out on scales much 
smaller than the anthropic range, and effectively we recover (3). Clearly, the only way to 
avoid this averaging effect is if rj > m p , and 

M 4 ~ (Ap D ) anth . (18) 

The last equation is to ensure that a significant range of values of 0/p is sampled in the 
anthropic range (ApD) an th ^ 10 3 po, so that changes in the slope of the potential are appre- 
ciable. Otherwise the distribution for V* will be almost flat. Thus, aside from the fact that 
the washboard potential is already a somewhat contrived example [19], Eq. (18) implies an 
otherwise unnecessary adjustment of the parameter M. 

In what follows, we shall only consider models where there is no such ad- hoc adjustment. 
In this sense, our predictions may not be completely inescapable, but they can be considered 
generic. The situation can be compared with the predictions of inflation that the density 
parameter is Q = 1 and the spectrum of density perturbations is nearly flat. It is certainly 
possible, in the context of inflation, to have an open universe with Q < 1, or to have 
a markedly non-flat spectrum of density perturbations. But to achieve this, additional 
parameters must be introduced and adjusted to the desired outcome. 



4 If M 4 3> H^m p rj 1 the slow roll condition is not satisfied today and the field <p will be in any one of the local 
minima of the washboard. With some generic requirements on the inflationary parameters, the minima 



will have equal a priori probability within the anthropic range [1]. 
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III. THE ANTHROPIC FACTOR 

We now consider the effect of the anthropic factor n civ in Eq. (2). The physical situation 
is rather different for positive and negative po, so we consider these two cases separately. 

For positive po, the main change introduced by n C i V is that the time of earliest galaxy 
formation tE G in the anthropic range (1) is effectively replaced by the time at which the 
bulk of galaxy formation occurs. This is because a few early birds will not make a difference 
once we apply the principle of mediocrity. More precisely, we should take into consideration 
that the morphology of some galaxies could make them less suitable for the development of 
civilizations, and therefore 

n civ {p D ) = J da n(a,p D ) N civ (a). (19) 

Here, a denotes the set of parameters characterizing the type of galaxy (e.g. its size, density, 
etc.), n(a,pi)) is the number density of such galaxies that form per co-moving volume in 
regions characterized by p^, and N civ (a) is the number of civilizations per galaxy of type a. 
Suppose that the above integral receives a dominant contribution from galaxies of type a G - 
Then 

n civ (p D ) oc n(a G ,p D ), (20) 

and the relevant time for anthropic considerations is the time at which this type of galaxies 
form, which we shall denote by t G . With the assumption of a flat prior V*, it was shown in 
[11, 39] that the most probable value for a positive p D is the one characterized by 

t D ~ t G . (21) 

This fact was used in order to explain the observed time coincidence 

t D ~ t - (22) 

The last relation follows from (21), assuming that stars and civilizations develop on a 
timescale not much greater than t G , and therefore t G is comparable to t , defined as the 
time when most civilizations make their first determination of po- 

Connected with the above discussion, there is a prediction of the anthropic approach, 
which can be checked by a combination of observations and theoretical analysis. In a not so 
distant future, our understanding of galactic evolution and morphology may improve to the 
point where we can tell with some confidence which galaxies are suitable for sustaining plan- 
etary systems similar to our own, where civilizations can develop. The anthropic approach 
to the cosmological constant problems (CCPs) predicts that the conditions for civilizations 
to emerge will be found mostly in galaxies that formed (or completed their formation) at a 
low redshift, z ~ 1. 

In the standard cold dark matter cosmology, galaxy formation is a hierarchical process, 
with smaller objects merging to form more and more massive ones. We know from observa- 
tions that some galaxies existed already at z — 5, and the theory predicts that some dwarf 
galaxies and dense central parts of giant galaxies could form as early as z = 10 or even 20. 
The fraction of matter bound in giant galaxies (M ~ 1O 12 M ) at z — 1 (~ 20%) is some- 
what less than that in objects of mass ~ 1O 9 M at z — 3, or in objects of mass ~ 1O 7 M 
at z = 5 [30]. If civilizations were as likely to form in early galaxies as in late ones, then 
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Eq. (21) would indicate that, for a typical observer, the cosmological constant should start 
dominating at a redshift zq > 5. The corresponding dark energy density, 

p D ~ (1 + z G ) 3 p , (23) 

would be far greater than observed. Clearly, the agreement becomes much better if we 
assume that the conditions for civilizations to emerge arise mainly in the types of galaxies 
which form at lower redshifts, zq ~ 1. 

We now point to some directions along which the choice of zq ~ 1 may be justified. One 
problem with dwarf galaxies is that if the mass of a galaxy is too small, then it cannot retain 
the heavy elements dispersed in supernova explosions. Numerical simulations suggest that 
the fraction of heavy elements retained is ~ 30% for a 1O 9 M galaxy and is negligible for 
much smaller galaxies [37]. The heavy elements are necessary for the formation of planets 
and of observers, and thus one has to require that the structure formation hierarchy should 
evolve up to mass scales ~ 10 9 M Q or higher prior to the dark energy domination. This gives 
the condition zq < 3, but falls short of explaining zq ~ 1. 

Another point to note is that smaller galaxies, formed at earlier times, have a higher 
density of matter. This may increase the danger of nearby supernova explosions and the rate 
of near encounters with stars, large molecular clouds, or dark matter clumps. Gravitational 
perturbations of planetary systems in such encounters could send a rain of comets from the 
Oort-type cloud towards the inner planets, causing mass extinctions 5 . 

Our own Galaxy has definitely passed the test for the evolution of intelligence, and the 
principle of mediocrity suggests that most observers may live in galaxies of this type. Our 
Milky Way is a giant spiral galaxy. The dense central parts of such galaxies were formed 
at a high redshift z > 5, but their discs were assembled at z ~ 1 or later [32]. Our Sun is 
located in the disc, at a distance ~ 8.5 kpc from the Galactic center 6 . If this situation is 
typical, then the relevant epoch to use in Eq. (23) is the epoch zq ~ 1 associated with the 
formation of discs of giant galaxies. 

The above remarks may or may not be on the right track, but we emphasize once again 
that if CCPs have an anthropic resolution, then, for one reason or another, the evolution 
of intelligent life should require conditions which are found mainly in giant galaxies, which 
completed their formation at zq ~ 1. 

In order to estimate n(ac,PD) in Eq. (20), we shall need a simple quantitative criterion 
to specify the relevant type of galaxies. The most important parameter characterizing a 



5 The cross-section for disruption of planetary orbits is much smaller, and it would take a rather substantial 
increase of the density for this process to become statistically important. A.V. is grateful to David Spergel 
for a discussion of this issue. 

6 It has been noted [31] that this distance is close to the corotation radius, where the orbital velocity of the 
stars coincides with the rotational velocity of the spiral pattern. In other words, the motion of the Sun 
relative to the spiral arms is rather slow, and as a result, the periods between spiral arm crossings are 
rather long (~ 10 8 yrs). Spiral arms are the primary sites of supernova explosions. They are also rich in 
giant molecular clouds, and are therefore very hazardous to life. It has been argued in [31] that spiral arm 
crossings are responsible for the major mass extinctions observed in the fossil record. Then one expects 
that habitable planetary systems are to be found mainly in the vicinity of the corotation radius, since 
mass extinctions at a rate much greater than once in 10 8 yrs may be too frequent for intelligent life to 
evolve. (Note that it took us 6.5 x 10 7 yrs to evolve since the last great extinction.) 
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galaxy is its mass M. For the Milky Way it is M MW ~ 10 12 M Q [33], and the above 
discussion suggests that we identify the relevant galaxies with gravitationally bound halos 
of this mass. (Note that this is also the typical mass of galaxies, which contain most of 
the luminous stars in the Universe.) It should be recognized, however, that the choice of 
this characteristic mass scale is somewhat uncertain, so we shall illustrate how our results 
are affected by choosing a larger or a smaller mass. 

Our Galaxy is a member of the Local Group cluster, whose mass has been estimated as 
[34] M LG ~ 4 x 10 12 Mq. It is conceivable that the gas captured in this cluster is later 
accreted onto the member galaxies and thus affects the properties of their discs. There 
seems to be no justification to consider larger mass objects, and we shall regard M LG as an 
upper bound on the potentially relevant mass scales. On the lower mass end, we shall use 
M ~ 10 11 M , which is roughly the mass of the bright part of our Galaxy, up to ~ 10 kpc 
from the center. (We note that M M w is probably a more reasonable choice, because the 
properties of the disc depend on the total mass of the halo [35].) 

We now consider negative po- The scale factor of a universe filled with nonrelativistic 
matter and dark energy with p D < is given by 

a(t) = sin 2 / 3 , (24) 

where to = (I/QttGIpdI) 1 ^ 2 ■ The matter density pu initially decreases while the universe 
expands, but at t — 7t£d/2, when it reaches the value pu = ~Pd, the universe stops its 
expansion and starts recontraction. The matter density grows in the contracting phase, and 
thus pm > \pd\ throughout the evolution. The structure formation in a universe with a 
negative po proceeds as usual until t ~ to, but then the growth of density perturbations 
accelerates during the contraction, so that all overdensities collapse to form bound objects 
prior to the big crunch. For to ^ to, giant galaxies will form at about the same time as 
they did in our part of the universe and will have similar properties (with a possible caveat 
indicated below). However, for t D <C t halos of the galactic size will be forced to collapse 
at a much earlier time t ~ to, and their density will therefore be much higher than that of 
our Galaxy. This would probably make such halos unsuitable for life. 

These considerations suggest that the anthropic factor effectively constrains to to be in 
the range 

to > to (25) 

for both positive and negative pn- There is, however, an additional factor that could make 
negative pa less probable. For p D > 0, structure formation effectively stops at t > tc, and 
the existing structures evolve more or less in isolation. This may account for the fact that 
discs of giant galaxies take their grand-design spiral form only relatively late, at z ~ 0.3. 
The discs are already in place at z ~ 1, but they have a very unsettled, irregular appearance 
[32]. On the other hand, for po < the clustering hierarchy only speeds up at t > to, and 
quiescent discs which may be necessary for the evolution of fragile creatures like ourselves 
may never be formed. 

Another factor to consider is the characteristic time £/ needed for intelligence to develop. 
For positive po, this factor is unimportant, since the time after the dark energy domination 
is practically unlimited, but for negative pr> the available time is bounded by t < nt D , and 
the effect of tj requires a closer examination. 

We first note that ti t is unlikely, since then it is not clear why it took so long for 
intelligence to develop on Earth. (The total time of biological evolution, from the origin of 
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life on Earth till present, is estimated at ~ 3.5 x 10 9 yrs.) For ij ^> to, we note that the main 
sequence lifetime of stars believed to be suitable to harbor life is t* ~ (5 — 20) x 10 9 yrs ~ to 
(see [11] for a discussion of this point). If tj ^> t ~ t+, most of these stars will explode as red 
giants before intelligence has a chance to develop. Carter [6] has argued that this is the most 
likely scenario 7 . In this case, the number N civ is suppressed by a factor ~ min{t*, to}/ti ~ 
to/ti, where we have used (25) in the last step. For positive pr>, the suppression is by a 
factor ~ t±/ti, which is of the same order of magnitude. 

We conclude that the precise value of tj has little effect on the relative probability of 
positive and negative pc,. If the accelerated clustering hierarchy is detrimental for life, then 
the probability for negative pu is suppressed; otherwise the two signs of pr, are equally likely. 
In either case, we should not be surprised that po is positive in our part of the universe. In 
the following sections we shall focus on the positive values of p D . 

IV. PREDICTION FOR THE EQUATION OF STATE 

A rather generic prediction of models where both CCP's are solved anthropically is that 
the equation of state of dark energy is given by po = wpo, with 

w = -l±l(T 5 . (26) 

The error bars correspond to the precision to which the observable universe can be aproxi- 
mated by a homogeneous and isotropic model. In models where px is the energy density of 
a four-form field, this equation of state is guaranteed by the fact that the four-form energy 
density is a constant and can only change by the nucleation of branes (other than that, 
it behaves exactly like an additional cosmological constant). If px is a generic scalar field 
potential, the slow roll conditions (6) are likely to be satisfied by excess, by many orders 
of magnitude, rather than marginally. For instance, for the quadratic potential (10), these 
conditions imply the constraint (11). It would be contrived to arrange for the condition to 
be satisfied marginally, since the whole point of the present approach is to have pa cancelled 
regardless of its precise value (which is not known to us even by order of magnitude). If 
the slow roll conditions are satisfied by excess by just more than three orders of magnitude, 
then the kinetic energy of the scalar field will be less than its potential energy by more than 
six orders of magnitude, and Eq. (26) follows. 

There are certainly models for dark energy, some of them with anthropic input, were 
(26) is not satisfied. For instance, Kallosh and Linde [13] recently considered a supergravity 
model where the time coincidence problem is solved anthropically, and where (26) does 
not hold. However, their model does not solve the old CCP, since it is assumed that the 
cosmological constant vanishes in the observable matter sector due to some unspecified 
mechanism. Likewise, (26) does not hold in the usual quintessence models [21], which have 
no anthropic input at all, but which do not address the CCP's [1, 19], or in models of 
k-essence [22], where only the time coincidence is partially addressed. 

A possibility worth discussing is the case of models where the slow roll parameters are 
themselves random fields. Consider, for instance, the following model: 

Pd =Pa + /i 2 (V# 2 . (27) 



7 The coincidence t] ~ is unlikely, since the evolution of life and evolution of stars are governed by 
completely different processes. 
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If the probability distribution for the new scalar field tp were such that all values of /i 2 are 
equiprobable, then one might imagine that the order of magnitude of \j? would be such that 
the slow roll conditions would be marginally satisfied. 8 



Currently favoured values for the dark energy density and for the Hubble parameter are 
Qn ~ .7 and h « .7 [24, 25], both with error bars of the order of 10%. While observations 
are not very accurate, we would like to challenge the status quo and boldly use the anthropic 
approach to the CCP's to make predictions for these two parameters. As we shall see, this 
approach predicts that Qd is likely to be somewhat higher, and that h is likely to be smaller 
than those currently favoured values. 

The basic reason why we expect Qn to be larger is the following [14, 15]. The growth 
of density fluctuations in a universe with a positive cosmological constant effectively stops 
at the redshift Zd when the cosmological constant starts dominating. This is given by 
(1 + z D ) ~ (VLo/VLm) 1 ^ , where Vl M = 1 — is the matter density parameter. According 
to (21), we expect zd ~ zq, where zq is the epoch when the relevant galaxies were formed. 
With zq ~ 1, this corresponds to (Q^/Qm) ~ 8, which in turn implies Qn ~ -9. (For 
zq > 1, we would obtain an even higher value for p^.) This prediction can be made more 



The new field ip must also be a light field and hence its prior distribution is in principle calculable. We 
can actually consider a more general form of the potential for tp and (p, 



Around any point (?po, (po) on the curve 7 defined by V(ip, (p) = —p\, the potential can be approximated 
by a linear function of the fields. Moreover, we can always rotate coordinates in field space so that tp is 
directed along 7, and cp is orthogonal to it, 



Here V$ is the gradient of the potential at that point. During inflation, both fields %p and <p are randomized 
by quantum fluctuations. Hence, the prior probability distribution is given by the area in field space 
T^dipdcp oc dtpdep, which leads to 



Along the curve 7, the values of V$ that will carry more weight per unit distance along the curve are 
those for which the slope is smaller. In the published version of this paper, it was (incorrectly) concluded 
from this observation that "given a model where the slope of the potential is variable, smaller values of 
the slope are preferred a priori, and there is no reason to expect that the slow roll conditions should be 
satisfied only marginally." However, this is not necessarily correct in general, since the slope could be large 
for a very large range of tp, and large slopes may end up dominating the integral in Eq. (29). A detailed 
analysis of this point has been given in [42] , where it is shown that there is a wide class of two-field models 
for which a small slope is favoured a priori (and for which the conclusions of Section IV do indeed apply) . 
However, there is an equally broad class of models for which a large slope is favoured a priori. In that 
case, the posterior probability distribution for the slope is such that the slow roll condition is satisfied 
only marginally, and some departure from the vacuum equation of state may be expected (see [43, 44]). 



V. PREDICTIONS FOR U D AND h. 



PD = PA + V(lp,(p). 



(28) 



PD ~ V0(V>O,0o)(</>- <Po)- 




(29) 
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quantitative [17] by using the distribution (2). As we shall see, the precise predictions 
depend not only on Qjj but also on h. 

Throughout this section, we shall assume that p D > as part of our prior. In a universe 
filled with pressureless matter and with a dark energy component p D > 0, the scale factor 
behaves as 

a(t) = sinh 2 / 3 ( — 

\tD, 

where to = (I/QttGpd) 1 ^ 2 - A primordial overdensity will eventually collapse, provided that 
its value at the time of recombination is larger than a certain value 5 r c ec . In the spherical 

collapse model, this is estimated as 8 r c ec (pD) = l.lSxlic, where x rec = x(t rec ) [27]. Here, we 
have introduced the variable 

The number of galaxies n(M,p£,) of mass M that will form per unit comoving volume in 
a region characterized by the value po of the dark energy density, is proportional to the 
fraction of matter that eventually clusters into this type of galaxies. In the Press-Schechter 
approximation [26, 27], this is given by 



n civ {Pn) oc n(M, p D ) oc erfc f^f^fL) . (31) 



V2a rec (M), 

Here, erfc is the complementary error function, and a rec (M) is the dispersion in the density 
contrast at the time of recombination t rec . As argued in the preceeding section, we shall 
assume that most civilizations are formed in galaxies characterized by a mass M ~ Mmw ~ 
1O 12 M (although we shall also consider slightly larger and smaller masses). 

The factor n civ depends on the parameter a rec , which in turn depends on the amplitude 
of density perturbations generated during inflation. The value of be inferred from 

the normalization of CMB anisotropies, but for this task, both the present value of Q D and 
the value of the Hubble parameter h would be needed. Since these are the parameters we 
wish to make predictions about, it would be somewhat contrived to use them at this point 
to make an inference about a rec . 

Another factor to consider is that a rec may be different in distant regions of the universe 
(where, as a consequence, galaxies would form earlier or later). In models where the inflaton 
field has only one component, the value of a rec is the same in all regions of the universe. 
However, if the inflaton field has more than one component, the amplitude of density per- 
turbations depends on the path followed by the inflaton on its way to the minimum of the 
potential. In such models, it is possible for a rec to vary over distances much larger than the 
presently observable universe. 

To make our discussion sufficiently general, we shall consider that a rec is itself a random 
variable with unspecified prior. This prior may be determined by processes occurring during 
inflation, or it may just reflect our ignorance of the actual value of the fixed parameter a rec . 
Then, Eq. (2) is generalized to 

dV(p D , a rec ) ~ n civ V*(p D , a rec )dp D da rec . (32) 

In this context, the generic expectation that the prior does not depend on p D in the anthropic 
range [see Eq. (3)], translates into 

V*(PD, CTree) ~ V*{o rec)- 
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dP/dln y 




Y 



FIG. 1: The distribution (34). 



Substituting (31) into (32), we have 

dV(p D , a rec ) oc erfc '- — V*(a rec ) dx rec da rec , (33) 

where we have used that f2^(t rec ) w 1 in all regions of interest, so that dpu oc dx rec . 
Introducing y = x rec a~^ c , the change of variables (x rec , a rec ) — > (y, cr rec ) produces a Jacobian 
proportional to of ec , and we have 9 dV(y,a rec ) ~ /(y)a 3 ec P*(a rec )ci?/<i<7 rec , where f(y) does 
not depend on a rec . Integrating over a rec leads to the normalized distribution 

dV(y) = (.80) 3 7r- 1/2 erfc(.80i/ 1/3 ) y d\ny, (34) 

which is uncorrected with a rec . 

The variable y can be expressed in terms of observable quantities, as we shall see below, 
and from (34) we should expect y ~ 1 by order of magnitude (See Fig. 1). More precisely, 



9 The appearance of the factor af ec in the posterior distribution was noted in Ref. [11]. This factor implies 
that the prior P*(cr rec ) should decay faster than a~^ c at large values of the density contrast. Otherwise, 
the posterior distribution is not normalizable (and we should expect both a large value of the effective 
vacuum energy and of the linear density contrast, in contradiction with observations). In the published 
version of the present paper this factor was omitted, due to an incorrect normalization of Eq. (33). We 
thank Takahiro Tanaka for drawing our attention to this point. 
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we expect y > .79 with probability 

P( y > .79) = .68 (la c.l), (35) 

and y > .07 with probability 

P( y > .07) = .95 (2cr c.l). (36) 

We shall denote these two equations as the la and 2a confidence level predictions for y. Let 
us now show how these translate into confidence level curves for the expected values of the 
parameters Qd and h. Here, and in what follows, Qd will denote the present value of the 
dark energy density parameter in our observable universe. 

Let us first express the "observed" value of y, which we shall denote as y , in terms of 
Q D and h. The density contrast at present is given by a = G(x ,x rec )a rec , where, assuming 
z rec ^> 1, the growth factor is given by [17] G(x ,x rec ) = x r lJ 3 F(Q D ), with 

F((l D ) = 5 -(l D ^ / *i (37) 

Therefore, 



Vo = • (38) 

The linearized density contrast at present a can be inferred from measurements of CMB 
temperature anisotropies, as described e.g. in [17, 28]. Since the spectrum is expressed as 
a function of wavelength, the mass scale has to be converted into a length-scale. A halo 
of mass M corresponds to a co-moving radius R(M) = (3M/47rp ) 1//3 - The mean matter 
density of the universe is given by p = 1.88 x 10~ 29 ^m^ 2 g/cm 3 , which leads to 

1/3 



R(M) = Mh-V^ Mpc. 



Assuming an adiabatic primordial spectrum of scalar density perturbations, characterized 
by a spectral index n, we have 

a (R) = {c im T)^l 2 5 H K l l\R). (39) 

Here, cioo = 2.9979 is the speed of light in units of 100km s _1 and 

T = n M hexp[-Q b (l + V2hVLl})} 

is the so-called shape parameter, with the density parameter in baryons. For numerical 
estimates, we shall take Vl h h 2 ks .02. The dimensionless amplitude of cosmological pertur- 
bations inferred from the COBE DMR experiment is given by [28, 29] 

5 H = 1.91 x iQ-5 exp[1.01(l-n)] ^ 80 -, 051n n M[1 _ lg(1 _ ^ _ (4Q) 
a/1 + r(.75 - .13Q 2 D ) 

The parameter r denotes the ratio of tensor to scalar amplitudes. Note that the effect of 
tensors is to make 5h a bit smaller (although not very significantly). Finally, 

POO 

K(R) = / q {n+2) T 2 (q)W 2 (q T hR M P c- l )dq, 
Jo 
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FIG. 2: Contours of the function yo(^Dj^) given in (38), corresponding to the la (lower curves) 
and 2cr (upper curves) predictions represented by Eqs. (35-36). The excluded region lies to the left 
of the curves. The thick solid lines assume that the dominant contribution to n c i v is in galaxies of 
mass M = Mmw = 10 12 M Q . For comparison, we show the predictions for different choices of the 
mass. The short dashed curves correspond to the mass of the local group Mlg = 4 x 1O 12 M , and 
the long dashed curves correspond to the mass of the bright inner part of our galaxy M = 1O 11 M0. 
A scale invariant spectrum of density perturbations is assumed. 

where the transfer function is given by T(q) = (2.34g) _1 ln(l + 2.34g)[l + 3.89g + (16. lq) 2 + 
(5.46g) 3 + (6.71g) 4 ] -1 / 4 and the window function is given by W(u) = 3w~ 3 (sin-u — ucosu). 

Substituting (39) in (38), and using Qd + Qm = 1, we obtain the function y = y (flr>, h). 
Contour lines of this function, corresponding to the la and 2a predictions represented by 
Eqs. (35-36), are plotted in Fig. 2, assuming that the dominant contribution to n C i V is in 
galaxies of mass M = Mmw = 10 12 M Q (thick solid lines). We also consider the predictions 
for different choices of the mass, as discussed in Section III. The short dashed curves cor- 
respond to the mass of the local group M LG = 4 x 1O 12 M , and the long dashed curves 
correspond to the mass of the bright inner part of our galaxy M = 1O 11 M . The effect of 
a tilt in the spectral index is plotted in Fig. 3. Both of these figures ignore the effect of 
tensor modes in the normalization (40). Tensor modes tend to lower the value of Sh, and 
hence they tend to make the bounds somewhat less stringent. The effect, however, is not 
dramatic. Even for r as large as .5, the effect on the curves is comparable to the effect of 
lowering the spectral index by .05. 

Expressions similar to Eqs. (34)- (38) were already contained in the exhaustive analysis of 
the problem given by Martel, Shapiro and Weinberg (MSW) in [17], where a rec was treated 
as a fixed parameter. However, our use of these expressions is somewhat different. MSW 
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h 




FIG. 3: Effect of a tilt in the spectral index of density perturbations. As in Fig. 2, the thick solid 
lines correspond to a scale invariant spectrum n = 1, and a mass M = Mmw = 1O 12 M0. The 
long dashed line and the short dashed lines correspond to tilted spectra, with n = .95 and n = .9 
respectively. 

noted that the existing observations indicate a value of £l D ~ 0.6 — 0.7 and used Eq. (33), 
with h = 0.7 to show that this range corresponds to probabilities from 2% to 12%, depending 
on the values chosen for the galactic scale M and the spectral index of perturbations n. They 
concluded that "anthropic considerations do fairly well as an explanation of a cosmological 
constant with [Qd] hi the range 0.6 — 0.7". However, one cannot help but feel disappointed 
by the somewhat low values of the probabilities. 

Our approach here is that anthropic models should be used as any other models - to 
make testable predictions. Thus, the goal is not so much to explain the value of Qd after it 
is determined by observations, but to predict that value at a specified confidence level. The 
contour lines in Figs. 2,3 indicate the la and 2a predictions of the model. If M ~ M MW 
proves to be the relevant mass giving the dominant contribution to n C i V , then the currently 
favoured model with Qd pa .7 and h ~ .7 is virtually excluded by the anthropic approach at 
the 2a level. Instead, this approach favours lower values of h and higher values of flo- 

These predictions can be turned around. If the values £l D < 0.7, h > 0.7 are confirmed 
by future measurements, then our model will be ruled out at a 95% confidence level, again 
assuming M ~ Mmw and a scale invariant spectrum. For a tilted spectrum, slightly lower 
values of are allowed at the same confidence level. The observational situation at the 
time of this writing is far from being clear. CMB and supernovae measurements yield [23, 24] 
Hd ~ 0.7 , while the observations of galaxy clustering give [38] Vt M = 0.18 ± 0.8, and thus 
^ D ~ 0.8. 
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VI. THE FUTURE OF THE UNIVERSE 

We finally discuss the anthropic prediction which is not likely to be tested any time soon. 
In all anthropic models, po can take both positive and negative values, so the observed 
positive dark energy will eventually start decreasing and will turn negative, and our part of 
the universe will recollapse to a big crunch. 

To be specific, we shall consider a scalar field model with a very flat potential. In the 
anthropic range (1), the potential can be approximated as a linear function, 

V{4>) » -V' <f>, (41) 

where V' is a constant and we have set = at V = 0. Once the dark energy dominates, 
the evolution is described by the usual slow roll equations 

3#0 = V'q, (42) 

# = ^VW, (43) 
where H = a/a and a(t) is the scale factor. The solution of (42), (43) is 

m = -041 - {t/u)f\ (44) 

a(t)=exp[47rm p - 2 (0 2 - 2 (t))], (45) 
where — 0* is the present value of and 

U = 87rt D (0,/m p ) 2 (46) 

is the time from the present to the beginning of recollapse. 

The slow roll condition (6) implies that 0* > m p . As we discussed in Sec. IV, we do not 
expect this condition to be only marginally satisfied, and thus 0* ^> m p . Then it follows from 
Eqs. (46) and (45) that t* ^> Sntp and therefore we should expect our region of the universe 
to undergo accelerated expansion for at least another trillion years before recollapse. 10 

The slow roll approximation breaks down at ~ — m p , so the above equations cannot be 
used to describe the evolution at > 0, where the potential becomes negative. A general 
analysis of models with negative potentials has been given in [36], where it is shown that at 
^> m p the dynamics becomes dominated by the kinetic energy of the field, 2 ^> |V(0)|. 
The corresponding evolution is described by 

</>(t) = —=L ln(t c - t) + const, (47) 

V67T 



a(t) oc (t c -t) 1/3 , (48) 

where t c is the time of the big crunch. The linear approximation (41) for the potential breaks 
down at sufficiently large 0, but in this regime the form of the potential is unimportant and 
Eqs. (47), (48) still apply. 



10 This is in contrast with the model of Kallosh and Linde [13] discussed in Section IV, where the universe 
is expected to recollapse within 10-20 billion years. 



17 



During the dark energy dominated expansion, the ordinary nonrelativistic matter is di- 
luted by the exponential factor (45). When the contraction starts, the density of matter 
begins to grow as pu oc (t c — t) -1 . However, the kinetic energy of the field grows much 
faster, 2 oc (t c — t)~ 2 , and thus ordinary matter forever remains a subdominant component 
of the universe. 

VII. CONCLUSIONS 

We now summarize the predictions that follow from the anthropic approach to the CCP's. 

(1) In the simplest models where the dark energy density takes discrete values, or where 
the dark energy density is due to the potential energy of a single scalar field, the dark energy 
equation of state is predicted to be that of the vacuum, 

Pd = wp D , (49) 

where w — — 1 with a very high accuracy. This distinguishes the anthropic models we 
discussed here from other approaches, such as quintessence [21] or /c-essence [22]. 11 

(2) The anthropic predictions for the dark energy density and for the Hubble parame- 
ter h are given in Figs. 2 and 3 of Section V 12 . We show the areas in the Vt d — h plane that are 
excluded at la and 2a confidence levels. The excluded areas depend on the assumed galactic 
mass M and on the spectral index n of the density fluctuations. For M = Mmw — 1O 12 M 
the currently popular values Vtu = 0.7, h = 0.7 are marginally excluded at 2a confidence 
level for a scale invariant spectrum n — 1. Lowering the spectral index relaxes the bounds 
somewhat. For h > 0.65 and n > .95, the la prediction is £Id > 0.79. These anthropic 
constraints get weaker when the relevant mass scale M is increased. For example, with 
M = 4 x 10 12 M Q a value as low as = 0.63 is still allowed at the 2a level for a scale 
invariant spectrum. The la prediction in this case is Qn > 0.78 (for h = 0.65). 

(3) Conditions for intelligent life to evolve are expected to arise mainly in giant galaxies 
that form (or complete their formation) at low redshifts, zq < 1. 

(4) The accelerated expansion will eventually stop and our part of the universe will 
recollapse, but, at least in the framework of the simplest models, it will take more than a 
trillion years for this to happen. Of course, this prediction is not likely to be tested anytime 
soon. 13 . 

The above predictions apply to models where both CCP's are solved anthropically. For 
comparison, we may consider other models. For instance, it is conceivable that a small value 
of the cosmological constant will eventually be explained within the fundamental theory. 
(We note the interesting recent proposal by Dvali, Gabadadze and Shifman [40] in this 
regard.) Even then, the coincidence problem will still have to be addressed. One possibility 



For a discussion of multifield models, see [42]. There is a class of two field models where the above 
prediction does not hold. These are the models where the prior distribution favours a large slope of the 
potential. In general, the equation of state pd = w(z)po depends on the redshift z, and increases with 
time as the dark energy fields pick up kinetic energy, on their way to negative values of the potential. The 
prediction for w(z) in the general case has been discussed in [43, 44]. 

These predictions are implicit in the earlier analysis by Martel, Shapiro and Weinberg [17] 

Again, this prediction is somewhat different in two field models where the prior distribution favours a 

large slope of the potential [42] 



18 



is that pd is truly a constant, while the amplitude of the density fluctuations a rec is a 
stochastic variable. With some mild assumptions about the prior probability distribution 
^♦(cv-ec)) ^ can be shown [1] that most galaxies are then formed at about the time of vacuum 
domination. In this class of models, predictions (1) and (3) still hold, while the other two 
predictions no longer apply. 

Another possibility has been recently discussed by Kallosh and Linde [13]. They assumed 
an M-theory inspired potential 

V{4>) = A(2 - coshV2(j)) (50) 

with a stochastic variable A. An interesting property of this potential is that its curvature 
is correlated with its height (at = 0). As a result, the universe tends to recollapse 
within a few Hubble times after the dark energy comes to dominate. Assuming that other 
contributions to the vacuum energy are somehow cancelled (that is, that the old CCP 
is solved by some unspecified mechanism), Kallosh and Linde argue that the coincidence 
to ~ ti is to be expected, where ti is the time it takes intelligent life to evolve (they assume 
it to be ~ 10 10 yrs). Predictions (l)-(3) are not applicable to this model. The model does 
predict recollapse of the universe, but the corresponding timescale (~ 10 10 yrs) is much 
shorter than the anthropic prediction (4). 

Anthropic arguments are sometimes perceived as handwaving, unpredictive and unfal- 
sifiable lore, of questionable scientific validity. In our view, the results presented in this 
paper should dispel this notion. Here, we have used the anthropic approach to make sev- 
eral quantitative predictions, some of which may soon be checked against observations. It 
should also be emphasized that, for the particular case of dark energy, there are at present 
no alternative theories explaining both CCP's, or making generic predictions of comparable 
accuracy. 

The present bound on the equation of state parameter w from the CMB and supernovae 
measurements is [41] w < —0.7, which is consistent with the anthropic prediction of w = — 1. 
The value of w = — 1 is usually associated with a plain cosmo logical constant. However, 
if in addition to this equation of state, observations confirm some of the other predictions 
presented above, this may be taken as an indication that the dark energy is dynamical. Thus, 
a better understanding of structure formation and galactic evolution may in fact reveal a 
crucial property of dark energy, with important implications for particle physics. 
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